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ABSTRACT 

Phase transitions induced by high temperatures and strong magnetic 
fields are investigated in the Standard model. The consistent effective po- 
tential including the one-loop and ring diagram contributions is calculated 
and investigated for the wide range of the fields, temperatures and Higgs 
boson mass m#. All other particles - fermions and bosons - are taken into 
account with actual values of their masses. This effective potential is real 
at sufficiently high temperatures . It is shown that symmetry restoration is 
a first type phase transition for m# < 70 Gev . For heavier Higgs particles 
and the field strengths H > 0.1 • 10 24 G the local electroweak minimum could 
not be realized at all. Hence, the upper limit on the Higgs mass as well as 
the limit on the field strengths in the phase transition epoch follow. 

1 .Introduction 

The concept of symmetry restoration at high temperature has been inten- 
sively used in studying the evolution of the universe at its early stages. Nowa- 
days it gives a possibility to investigate various problems of cosmology and 
particle physics [JJ,0]. In particular, the type of the electroweak phase tran- 
sition and hence a further evolution of the universe depends on the mass m# 
of the Higgs boson. Most investigations of the electroweak phase transition 
have included into consideration a high temperature environment as the main 
ingredient 0,0. But in recent a few years cogent arguments followed from 
different approaches in favour of the presence of strong magnetic fields at 
that stage have appeared M , |6| ( recent survey on the magnetic fields in the 
universe is Ref.0). So, the phase transition at high temperature and strong 
fields has to be of interest . Moreover, at present time when all masses of 
fundamental particles, except m#, are known it is possible to investigate in 
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details the phase transition as the function of this parameter and to deter- 
mine the properties of the vacuum, its structure and to specify a so-called 
metastability bound on m H 

One of the ways to have strong magnetic fields in the electroweak phase 
transition epoch was proposed by Vachaspati ||. From his analysis it fol- 
lows that under very general conditions the fields H ~ Tf in the patches of 
sub-horizon scales can be generated during a large class of grand unified tran- 
sitions where Tj is transition temperature. The second one is formation 
of the Savvidy vacuum magnetic state at high temperature (H ~ gT 2 , g is 
gauge coupling constant) [^3|,||7|,[pl[. In latter case only the abelian field 
configurations could arise spontaneously since they are sourceless. For many 
problems it is important to estimate the field strengths presented, but it is 
difficult to realize that without detailed investigations within specific models. 
Usually, only one type of fields is considered. Therefore, results obtained in 
such a way give an upper estimate of the field. This remark is relevant to our 
present investigation. In current literature along with usual magnetic field 
a hypercharge magnetic field presenting in the restored phase is discussed ( 
see for example 0). The latter one is converted into ordinary magnetic field 
during the phase transition. 

In what follows, we will consider the case when the magnetic field is 
presenting in both broken and restored phases. We believe that this is a good 
approximation which gives possibility to investigate essential features of the 
phase transition. This scenario assumes that the field has been generated 
at a GUT scale via the Savvidy mechanism and presented during the phase 
transition. Such a picture is likely since, as it follows from our consistent 
calculation of the effective potential which found to be real at sufficiently high 
temperatures ( and, in particular, in the restored phase), constant magnetic 
field is stable. This is the most important point of the present analysis. In 
any case, it may have relevance to the description of the phase with broken 
symmetry at high temperatures and strong fields. 

The discussed "primordial" fields are usually considered as seed fields 
responsible for generation of the observed magnetic fields in galaxies J7J. 

Various aspects of the phase transitions in magnetic fields at high tem- 
perature have been investigated by many authors H- p5|| . In Refs. [12], [19 



the phase transitions derived with the effective potential (EP) of the bosonic 
part of the Salam- Weinberg model and the vacuum structures of the phases 
have also been described. In Ref.fT5|| in addition to the temperature and 
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magnetic fields a chemical potential was incorporated. But the role of the 
fermions has not been investigated in detail. However, due to a rather heavy 
i-quark mass, m t ~ 175 Gev, an unbounded global minimum of the EP is 
produced in addition to the electroweak local one for not very heavy Higgs 
scalars ||. As will be shown, strong magnetic fields influence essentially 
the phase transition dynamics realized in this case. Another aspect of the 
electroweak phase transition, which also was not investigated but plays an 
important role, is the influence of so-called ring diagrams at high tempera- 



ture and strong field. At zero field it was investigated in Refs. ||16| J171| where 



it has been shown the importance of these diagrams for determining the type 
of the phase transition. In Ref.|l7| the t-quark mass was chosen of order 
llOGev. So, taking into account the present day data, it should be consid- 
ered as a qualitative estimate of the role of ring diagrams even for zero-field 
case. 

In the present paper we study the electroweak phase transition at high 
temperatures and constant strong magnetic fields H . We calculate and 
investigate the one-loop effective potential(EP) and the contributions of ring 
diagrams. In contrast to previous considerations we include the content 
of all bosons and fermions with the corresponding masses. So, the only 
free parameter remains the mass m#. The role of ring diagrams is of great 
importance because they contain the terms which cancel out the imaginary 
part of the one-loop EP. The total EP is real at sufficiently high temperatures 
and suitable for investigations of symmetry behaviour. As it will be shown, 
the influence of fermions (mainly heavy quarks) is very essential. For values of 
mass m H < 70Gev and weak field strengths the electroweak phase transition 
is of the first type one. It is turned out that the magnetic field stimulates 
either the generation of the electroweak minimum and tends to remove the 
potential barrier separating this metastable state from the unbounded global 
minimum produced due to heavy fermions. At any high temperature there 
exists the corresponding field strength H at which the metastable vacuum 
could not be produced at all. This property can be used to find an upper limit 
on H at the transition epoch. The paper is organized as follows. In Sects. 2, 3 
the one-loop contributions of bosons and fermions to the EPV^(T, H, C ) are 
calculated in the form convenient for numerical study. In Sect. 4 we compute 
the contributions of ring diagrams. Further in Sect. 5 symmetry behaviour 
is investigated for a number of values of m# and H . For comparison, we 
consider separately the cases when only the one-loop EP is taken into account 
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and when the ring diagrams are included. Discussion of the results obtained 
is given in Sect. 6. 



2. Boson contributions to V^(T : H : (j} c ) 

The Lagrangian of the boson sector of the Salam- Weinberg model is 

+T ($+$) "I ($+$)2 ' (1) 

where the standard notations are introduced 

F; v = d,A:-d l/ A; + ge abc A b IM A c l/} 
G^ u — dpB v — d u B^, 

D, = d, + UgA^ + Ug'B,. (2) 
The vacuum expectation value of the field $ is 

The fields corresponding to the W^-,Z-bosons and photons, respectively, are 
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A » = ^=v(9'Al + gB»). (4) 

The external electromagnetic field is introduced by splitting the potential 
in two parts: A^ = A^ + A R , where A R describes a radiation field and 
A = (0, 0, Hx 1 , 0) corresponds to the constant magnetic field directed along 
the third axis. We make use of the gauge- fixing conditions 



2 / 2 

d,W^ ± leA^ =F = C^x), (5) 
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d^- l -{g 2 + g' 2 ) 1/2 <p z = C z , (6) 

where e = gsin9 w ,tang9 w = g' / g, 4> ± , 4> z are the Goldstone fields, £, £' are 
the gauge fixing parameters, C , C z are arbitrary functions and <p c is a scalar 
condensate value. In what follows, all calculations will be done in the gen- 
eral relativistic renormalizable gauge (||),(@) and after that we set = 
choosing the unitary gauge. 

To compute the EP in the background magnetic field let us introduce 
the proper time,s, representation for the Green functions 

oo 

G ab = -if dsexp{-tsG- iab ) (7) 
o 



and use the method of Ref. [18], allowing in a natural way to incorporate the 



temperature into this formalism. A basic formula of Ref.|L8| connecting the 
Matsubara Green functions with the Green functions at zero temperature is 
needed, 

+oo 

Gf{x,x'-T) = ^(-l^+^Gf (x - [x]f3u,x' - n/3u), (8) 

— oo 

where Gjfis the corresponding function at T — 0, /3 = 1/T,u = (0,0,0,1), 
the symbol [x] means an integer part of x^//3, = 1 in the case of physical 
fermions and = for boson and ghost fields. The Green functions in the 
right-hand side of formula (§) are the matrix elements of the operators Gk 
computed in the states | x', a) at T = 0, and in the left-hand side the oper- 
ators are averaged in the states with T ^ 0. The corresponding functional 
spaces U° and U T are different but in the limit of TB0 U T transforms into 
U°. 

The one-loop contribution into EP is given by the expression 

= ^Tr\ogG ab , (9) 

where G ab stands for the propagators of all the quantum fields W ± , ± , ... in 
the background magnetic field H. In the s- representation the calculation of 
the trace can be done in accordance with formula ETI 



30 ds 



/ds 
—trexpi-isG^ 1 ) (10) 
s 
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Details of calculations based on the s-representation and the formula @ can 
be found, for example, in Refs.[lS],[19|,[24|. The terms with n = in Eqs.@, 
(|S[) give zero temperature expressions for Green's functions and effective po- 
tential V^ 1 ', respectively. They are the only terms possessing divergences. 
To eliminate them and uniquely fix the potential we use the following renor- 
malization conditions at H, T = 0|19]: 



d 2 V{<p,H) 
dlP 

dcj) 
d 2 V{(j),H) 



\H=o,<f>=s(oy- 



1 

2' 



\H=0,, 



=5(0)= 



lh=o,. 



=5(0) : 



m 



(11) 
(12) 
(13) 



where V{<j>, H) = V® + V™ + • • • is the expansion in a number of loops and 
5(0) is the vacuum value of a scalar field determined in a tree approximation. 

It is convenient for what follows to introduce the dimensionless quantities: 
h = H/H (H = M 2 w /e)A = <f>c/5(0),K = m 2 H /Ml B = (3M w ,r = 1/B = 
T/M W ,V = VI E 2 and M w = ^5(0). 

After reparametrization the scalar field potential is directly expressed in 
terms of the ratio K, 



y(°) 



}i 2 



+ K{- 



+ 



(14) 



2 v 4 

The renormalized one-loop EP is given by the sum of the functions 

V 1 = V (0) +V {1) ((j),H,K) + J 1 \ ( J),h,K,T), (15) 

where is the one-loop EP at T = 0, which has been studied already in 
Ref . ||23|| . It has the form: 



v« = v« + v 



(1) 



(16) 



where 
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- ^h 2 log(2h)} 

+ —{-2h 2 + (h 2 + h<p 2 )log(h + 2 ) + (h 2 - h(j) 2 )log \ h - 2 |] 

27T 

+ i^ah(<p 2 - h)9(h- (p 2 ), (17) 



(!) _ ^ „„-„2/j (±1 , 1 i4 
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aK 2 u 9 A , 3 i2 1., 30 2 -l, 27 j4 21,.. ,, . 

and u/ 1 ) is the temperature dependent contribution to the EP determined by 
the terms of formulae ©,© with n^O. 

We outline the used calculation procedure considering the 1^-boson con- 
tribution as an example ||24|1 , 



«> = ^ 7 % e-^)[ 1 + 2c ° S2g l ±eMihB 2 n 2 /4s). (19) 
J s 1 L sin s J i 

o 1 



111 



1 + 2 cos 2s 



As Eq.(|l7]), this expression contains an imaginary part for h > <fi 2 appeared 
due to the tachyonic mode e 2 = p\ + M 2 — eH in the H^-boson spectrum . 
It can be explicitly calculated by means of an analytic continuation taking 
into account the shift sQs — iO in the s- plane. Fixing <ft 2 /h > 1 one can 
rotate clockwise the integration contour in the s-plane and direct it along 
the negative imaginary axis. Then, using the identity 

1 oo 

_ L_ = 2 y e - s(2p+1) (20) 

sinhs ^ 

and integrating over s in accordance with the standard formula 

b 



J dss 11 - 1 exp(-- - as) = 2(^) n l 2 K n (2^d>), (21) 
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b > 0, one can represent the expression ([19]) in the form 

ReuW = -4-A(3^o w 2 ), (22) 

7T iD 



where 



x x> x 



u = ]T J2 — R i (nBx p ) ;x p = ((f) 2 + h + 2ph) l ' 2 (23) 



p=0n=l 71 



OC 



Wl = £ -K 1 (nBy),y= (0 2 - fc) 1/2 ( 24 ) 

n=l 71 

and in the region of parameters (f) 2 < h after analytic continuation 

7T °° I V I 

^-^E^iMUl), (25) 

Z n=l n 



^ 2 = Y, -K 1 {nBz),z = {<P 2 + h) 1 ' 2 , (26) 



and /^(x), Y n (x) are the Macdonald and Bessel functions, respectively. The 
imaginary part of u$ is given by the expression 

Imu 2 = -2a- £ ULlj^nB \ y |), (27) 
J n (x)is Bessel function. As is well known, the imaginary part of EP is sig- 



naling the instability of a system plfl. In what follows we shall consider 



mainly the symmetry behaviour described by the real part of the EP . As 
the imaginary part is concern, it will be cancelled in consistent calculation 
including the one-loop and ring diagram contributions to the EP. 

Carrying out similar calculations for the Z- and Higgs bosons, we obtain 

& 



a ™ 0" ,n 



7r ^— J cos w n 2 B 2 cos 

&Y 2 (nB\t\))- < 29 » 

n=l 
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where the variable t = [K w ( — 2 ~ 1 ' > )] 1 ^ 2 at 30 2 > 1 and series with the function 
Y 2 (x) has to be calculated at 30 2 < 1. The corresponding imaginary part is 
also cancelled as it will be shown below. 

The above expressions (|i6|) ,(p2|),(^8D,(^9|) will be used in numerical study- 



ing of the symmetry behaviour at different H, T. There is a cancellation of a 
number of terms from the zero-temperature contributions given Eqs. flT6|) and 
T-depended ones. This fact has a general character and was used in checking 
of the correctness of calculations. 

3. Fermion contributions to V^\H : T : (j} c ) 

To find the explicit form of the fermion contribution to the EP let us consider 



the standard unrenormalized expression written in the s -representation |22 
VP = --f- V (-If / %e^ m2s+p2n2 / As \eHs)cotheHs, (30) 

87r n=-oo Si 

m is a fermion mass. Here, we have incorporated the equation (^) to include 
a temperature dependence. In what follows, we shall take into account the 
contributions of all fermions - leptons and quarks - with their masses known 
at present time. Usually, considering a symmetry behaviour without field 
one restricts himself by a t-quark contribution only. But in the case of an 
external field applied this is not a good idea, since the dependence of on 
if is a complicated function of the ratio m'j/eH. So, at some fixed values of 
H, T different dependencies on H will contribute for fermions with different 
masses. Hence, a very complicate dependence on the field takes place in 
general. We include this in a total, carrying out a numerical calculations 
and summing up over all the fermions. Now, separating a zero temperature 

+0O oo 

contribution by means of the relation — 1 + 2 an d introducing the 

— oo 1 

parameter Kf = m'j/M^ = r ^G\ ukawa j ' g 2 we obtain for the zero temperature 
fermion contribution to the dimensionless EP , 

V/(M)= t E^(-20 2 + |^ 4 - 
/ z 

a V^/ 2 h2 i 2 I If 

- -E(?/ir lo s 



7T y y ^ J 6 ° K f 
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^ [29 > 2l ° gri( 2^ ) + (2C ' ( - 1} " 



6 

2 | 9/ | ft 



(31) 



where g/ is a fermion electric charge, the sum = (leptons)+3 (quarks) 

f /=i /=i 
counts the contributions of leptons and quarks with their electric charges. 

The Ti function is defined as follows (see, for example, survey p3|): 



1 1 1 

logTi(x) = J dylogT(y) + -x(x - 1) - -xlog(27r). 



(32) 



The finite temperature part can be calculated in a way, described in the 
previous section. In the dimensionless variables it looks as follows: 



-/ = -4 E{EEH 

H f p=0n=l 



(2ph + Kf^f^h 



Bn 



Kxifrph + K f 4> 2 ) 1/2 Bn) 
{^m + K^ + 2)h + Kf ^yn Bn)] } (33) 



Again, a number of terms from Eqs. ( plj) and (|33| ) are cancelled being summed 
up, as in the bosonic sector. 

These two expressions and the boson contributions obtained in Sect. 2 will 
be used in numerical investigations of symmetry behaviour. More precise, we 
consider the difference V'(H, T, 0) = ReV(H, T, 0) -ReV(H, T, = 0) giving 
possibility to determine the symmetry restoration. We will investigate the 
EP of two types - the one-loop contribution and the sum of that and the ring 
diagrams which are the next to leading order corrections at high temperature. 



4. Contribution of ring diagrams 

It was shown by Carrington [[17] that at T ^ a consistent calculation of the 
EP based on generalized propagators, which include the polarization opera- 
tor insertions, requires that ring diagrams have to be added simultaneously 
with the one-loop contributions. These diagrams essentially affect the phase 
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transition at high temperature and zero field flIE|,||17||. Their importance at 



T and H ^ was also pointed out in literature 
know, this part of the EP has not been calculated, yet. 
As is known 



14], [15] but, as far as we 
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the sum of ring diagrams describes a dominant contri- 
bution of great distances. It differs from zero only in the case when massless 
states appear in a system. So, this type of diagrams has to be calculated when 
a symmetry restoration is investigated. To find the correction V r i ng (H,T) at 
high temperature and magnetic field the polarization operators of the Higgs 
particle, photon and Z-boson at the considered background have to be calcu- 



lated. Just these calculations have been announced in Refs.[[L4]], fl5||. Then, 



V r i n g(H,T) is given by series depicted in figures 1,2. 



Figure f: Scalar fielck ring ^liagrams x gTving ^contribution to the effective po- 



tential / + i j 



Figure'z: Photon ana ^boson ring diagrams giving contribution to the 
effective potential" C J + C. J + ' ' ' 

Here, a dashed line describes the Higgs particles, the wavy lines represent 
photons and Z-bosons, the blobs represent the polarization operators in the 
limit of zero momenta. As also is known [T7) , in order to calculate the contri- 
bution of ring diagrams not the complete polarization operators U lll/ (k, T, H) 
but only their limiting expressions at zero momenta, IIoo(fc = 0,T,H), are 
sufficient. This limit, named the Debye mass, can be calculated from the EP 
of the special type. This fact considerably simplifies our task. 

Now, let us turn to calculations of V r i ng (H, T). It is given by the standard 



11 



expression |i~6f , fT7|j , [|T 



V ring = --^-Tr{[M 2 (0) + IMO)] 3 / 2 - M 3 (0)}, (34) 

where trace means the summation over all the contributing states, M(0) is a 
tree mass of corresponding state and n o(0) = U(k = 0,T,H) for the Higgs 
particle and n o(0) = Hoo(k = 0,T, H) are the zero-zero components of the 
polarization operators in the magnetic field taken at zero momenta. The 
above contribution has order ~ <7 3 (A 3 ) in coupling constant whereas the two- 
loop terms are to be of order ~ g 4 . As Hoo(0) the high temperature limits 
of polarization functions have to be substituted which have the order ~ T 2 
for leading terms and ~ g^cT, (gH) 1 '^, C , (gH) 1 ^ 2 /T « 1 for subleading 
ones. 

For the next step of calculation, we remind that the effective potential 
is the generating functional of the one-particle irreducible Green functions 
at zero momenta transferred. So, to have 11(0) we can just calculate the 
second derivative with respect to <p of the potential V^'(H, T, <p) in the limit 
of high temperature ,T >> 0, (eif) 1//2 and then set = 0. This limit can be 
calculated by means of the Mellin transformation technique (see for example 
[ P4| ) and the result looks as follows: 

6 cos A 
a ( 3A0g-^ 2 (O) )3/2T] 

3 2 

- l// c+g Hr-T + leHT{\^ + eH)^ 

+ l-eHT^tf - eHY>\ (35) 

3 3 

The parameter Ct — J2 G\ + 3 J2 Gf determines the fermion contribution 

i=i i=i H 

of the order ~ T 2 having relevance to our problem. We also omitted ~ T 4 
contributions to the EP. The terms of the type ~ log[T/ f((p, H)] cancel the 
logarithmic terms in the temperature independent contributions flT5|) , fl3T)|) ■ 
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Considering the high temperature limit we restrict ourselves by linear and 
quadratic in T terms, only. 

Now, differentiating this expression twice with respect to and setting 
6 = 0, we obtain 



n#(o) 



d 2 V^((j),H,T) 

W 10=0 

) 1/2 r,., 1 1. 



/ 



+ , (eH L ^ 2 (3V2C(-^)-l). (36) 
Sir sin y^p 2 2 



The terms ~ T 2 in Eq. ( |3~6"D give standard contributions to temperature mass 
squared coming from the boson and fermion sectors. The if -dependent term 
is negative since the difference in the brackets is 3v^2~C( — \i §) ~~ 1 — — 0; 39. 
Formally, this may result in the negativeness of the n(O)0 for very strong 
fields (eif) 1//2 > T. But this happens in the range of parameters where 
asymptotic expansion is not valid. Substituting expression (|36"D into Eq.(|3~|) 
we obtain (in the dimensionless variables), 



vU = —{(^k + n,(0)} 3/! + y<(^r-- (37) 



As is seen, the last term of this expression cancels the fourth term in Eq. (|35| ) , 
which becomes imaginary at 3<j) 2 < 1. This is the important cancellation 
preventing the infrared instability at high temperature. 

Before to proceed, let us note that Eq.fl3"5|) contains other term (the last 
one) which becomes imaginary for strong magnetic fields or small 2 . It 
reflects the known instability in the M^-boson spectrum which is discussed for 



many years in literature (see papers [19],[14|,|15], [24| and references therein). 



But it also will be cancelled out when the contribution of ring diagrams with 
the unstable mode is added. 

To find the if-dependent Debye masses of photons and Z-bosons the 
following procedure will be used. We calculate the one-loop EP of the W- 
bosons and fermions in a magnetic field and some "chemical potential",^, 
which plays the role of an auxiliary parameter. Then, by differentiating 
them twice with respect to fi and setting fi = the mass squared m 2 D will 
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be obtained. Let us demonstrate that in more detail for the case of fermion 
contributions where the result is known. 

The temperature dependent part of the one-loop EP of constant mag- 
netic field and a non-zero chemical potential induced by an electron-positron 



vacuum polarization is 22 



1 oo ™ J —R 2 1 2 

V« =- V(-l) m / ^exp(^ m 2 s)eHscoth(eHs)cosh{pifi), 

1 ' 4tt z ,i J s A 4s 

(38) 

where m is the electron mass, e = gsin8 w is electric charge and a proper-time 
representation is used. Its second derivative with respect to [i taken at /i = 
can be written in the form, 



«2ta(1) tt o oo °9 j 

= ^^(-l)^ 1 / -exp(-m 2 s -(3 2 l 2 /As)coth{eHs). (39) 

Expanding coth(eHs) in series and integrating over s in accordance with 
formula ( pij ) we obtain in the limit of T >> m, (eif) 1 / 2 : 

fX-ir^K^l) + l^l Kx {mpt) + ■■■} (40) 

Series in I can easily be calculated by means of the Mellin transformation 
(see Refs.|pH,P!5|). To have the Debye mass squared it is necessary to dif- 
ferentiate Eq . (p9|) with respect to [3 2 and to take into account the relation of 
the parameter [i with the zero component of the electromagnetic potential : 



jjQieAo |L4]. In this way we obtain finally, 



m 2 D = g 2 szn 2 9 w (Y ~ ^m 2 + 0((m/3) 2 , {eH(3 2 ))). (41) 

This is the well known result calculated from the photon polarization operator 
ZIJ. As one can see, the dependence on H appears in the order ~ T~ 2 . To 



find the total fermion contribution to m 2 D one should sum up expression (|4l|) 
over all fermions and substitute their electric charges. 

To calculate m 2 D for Z-bosons it is sufficient to take into account the 
fermion couplings with Z- field. It can be done by substituting fiQi(g/2cos9 w + 
gsin 2 9 w ) and the result differs from Eq. ( f41"D by the coefficient at the 
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bracket in the right-hand side which have to be replaced, g 2 sin 2 9$g 2 ( 4co ^ id + 
tang 2 9 w ). One also should add the terms coming due to the neutral currents 
and the part of fermion-Z-boson interaction which is not reproduced by the 
above substitution: 

m D = F^H 1 + 4smX)T 2 . (42) 
8cos z w 



Now, let us apply this procedure for the case of the J¥-boson contribution. 
The corresponding EP (temperature dependent part) calculated at non-zero 
T, fj, in the unitary gauge looks as follows, 

VP = HhiE / ^exp(-m 2 s-l 2 p 2 /4s)[ , * - +Asinh{eHs)]cosh{(3ln). 
o7r J s stnh{eHs) 

(43) 

All the notations are obvious. The first term in the squared brackets gives the 
triple contribution of the charged scalar field and the second one is due to the 
interaction with a W^-boson magnetic moment. Again, after differentiation 
twice with respect to /i and setting /j, — it can be written as 

d 2 V^ eH o2 d ™°fds . m 2 s l 2 [3 2 eH„. 3 A . u 

-at = w^wlkl T e ^ ( -^-^ )[ ^(^) +4sm/i(eifs)] - 

(44) 

Expanding sink x s in series over Bernoulli's polynomials, 

— E^i-^, (45) 



sinks s , n k\ 

and carrying out all the calculations described above, we obtain for the W- 
boson contribution to m 2 D of the electromagnetic field, 

m 2 D = 3g 2 sin 2 e w []-T 2 - ^-T(m 2 + gsin6 w H) 1/2 - J—(gsin6 w H) 

O All 07T 

+0{m 2 /T 2 ,{gstn6 w H/T 2 ) 2 )]. (46) 

Hence it follows that spin does not contribute to the Debye mass in the 
leading order. Other interesting point is that the next to leading terms are 
negative. 
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The contribution of the VT-boson sector to the Z-boson mass m? D is given 
by expression ( [Hf) with the replacement g 2 sin 2 9 w Qg 2 cos 2 9 w . 

Summing up the expressions (^Tj) and ( Pf ) and substituting them in 
Eq. (|34]) , we obtain the photonic part V^ ing where it is necessary to express 
masses in terms of the vacuum value of the scalar condensate C . In the same 
way the contribution of Z-bosons V* ing can be calculated. The only difference 
is an additional mass term of Z -field and the additional term in the Debye 
mass due to the neutral current ~ u^^uZmu . These three fields - 4>,J,Z 
,- which becomes massless in the restored phase, contribute into V r i ng (H,T) 
in the presence of the magnetic field. At zero field there is also a term due 
to the VT-boson loops in Figs. 1,2 . But when H ^ the charged particles 
acquire ~ eH masses. The corresponding fields remain short-range ones in 
the restored phase of the vacuum and therefore do not contribute. 

A separate consideration should be spared to the tachyonic (unstable) 
mode in the VT-boson spectrum: = p\ + M 2 — eH. First of all we note 
that this mode is produced due to a spin interaction and it does not influence 
the Goo(k) component of the P^-boson propagator. Secondly, in the fields 
eH ~ the mode becomes a long range state. Therefore it should be 
included in V r i ng (H,T) side by side with other considered neutral fields. But 
in this case it is impossible to take advantage of formula ([34]) and one has to 
return to the initial EP with generalized propagators . 

For our purpose it will be convenient to use the expression for the general- 
ized EP written as a sum over the modes in external magnetic field fl4jl, [|l5 



v& = || £ / ^ £ log[ ^ + + n(T ' Hm (47) 



271(3 t 



n=0,o-=0,±l 



where u t = , e 2 n = p§ + M 2 + (2n + 1 - 2a) eH and U(H, T) is the Debye 
mass of VT-bosons in a magnetic field. Denoting as Dq 1 (ps, H.T) the sum 
u)f + e 2 , one can rewrite eq. ( f47|) as follows: 



2*p ^ l 2* ; 

TT +00 +°° J 

+ / ^{log[l + (u;f+ P l + M 2 -eHr^(H,T)] 

°°-infty 
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+ £ log[l + D (el,H,T)U(H,T)}}. 



(48) 



Here, the first term is just the one-loop contribution of W-bosons, the second 
one gives the sum of ring diagrams of the unstable mode ( as it can easily be 
verified by expanding the logarithm into a series ).The last term describes 
the sum of the short range modes and should be omitted. 

Thus, to find V£-£J taWe one has to calculate the second term in Eq. fl48|) . 
In the high temperature limit we obtain: 

pR 

V rZi Me = - eH + m,T)) 1/2 - {Ml - eH)^}. (49) 

By summing up the one-loop EP and all the terms V r i ng , we arrive at the 
total consistent in leading order effective potential. 

Let us note the most important features of the above expression. It is 
seen that the last term in Eq.fl4*9|) exactly cancels the dangerous term in 
Eq . (|35|) . So, no instabilities appear at sufficiently high temperatures when 
U(H,T) > Ml — eH and the EP is real. To make a quantitative estimate 
of the range of validity of the total EP it is necessary to calculate the mass 
operator of iy-boson in a magnetic field at finite temperature and hence to 
find Il(if , T) . This is a separate and enough complicated problem which will 
not be solved here completely. Instead that below we restrict ourselves by 
the contribution to II (T) of the neutral Higgs particles only which can easily 
be calculated to give U(T) Higgs = j^g 2 T 2 . Since other particles have also to 
contribute into the temperature mass in leading T 2 order, the obtained mass 
to be lower value which can be substituted into V™™ g table . Just this value will 
be used in the following estimations. 



5. Symmetry behaviour in a magnetic field at high tem- 
perature 



Having obtained the one- loop EP described by formulae flTED , (P2|) , (P5|)-(P^) 
and the ring diagram contributions V r i ng we are going to investigate symmetry 
behaviour at high temperature and strong magnetic fields. We shall present 
the results in two steps. First, we consider the sum of the tree and one- 
loop effective potentials as the function of <p 2 at various fixed H, T and K 
. Then, we shall add the term V r i ng and calculate symmetry behaviour for 
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the total EP at the same fixed H, T, K. This will help to clarify the role 
of the plasmon contributions. Since constructed EP includes as an input all 
fundamental particles, we shall obtain new information about the electroweak 
phase transition. The one-loop EP contains the imaginary parts in some 
domains of <fi 2 . But these terms occur to be cancelled by the corresponding 
ones from V ring . Thus, only the real part, .Re V^ 1 ), is of interest. 

As usually fl23j , to investigate symmetry behaviour let us consider the 
difference V = Re[V(h, <p, K, B) - V{h,<f> = 0,B)] which gives information 
about the symmetry restoration. Below, we consider the case when the mass 
mil is equal to M w . Typical curves for small fields h and different values of 
B are plotted in Fig. 3. 

It is seen that the well known symmetry restoration (for the heavy fermion 
case) takes place. There are two minima - local, produces due to the Higgs 
mechanism, and and global one, generated by heavy fermion contributions. 
At low temperatures (big B) , the local metastable minimum is disposed near 
the value <fi 2 = 1 that corresponds to the spontaneously broken symmetry. 
With a temperature increasing the local minimum becomes shallower and at 
B ~ 0.1 removes to the value 2 = that signals the symmetry restoration. 
We see that typical second type phase transition takes place for K = 1. At 
the same time, the barrier separating two minima is increasing in height and 
width and so a tunneling to the global unbounded from below minimum is 
suppressed. 

In Fig. 4 we present the influence of the field on the symmetry behaviour 
at low temperature. As is seen, an increase in h leads to the getting deeper 
of the local minimum and to the growing up the barrier which separates 
the minima. In this way the magnetic field prevents tunneling to the global 
unbounded minimum. 

Now, let us investigate symmetry behaviour at high temperature and 
strong magnetic fields. The result of calculations is shown in Fig. 5. From 
the plot it follows that the field tends to decrease the temperature and stim- 
ulates the symmetry breaking. First the metastable vacuum is generated 
with the field increasing. As is seen, this is a homogeneous transition. When 
the strength h is growing further the potential barrier separating the local 
and global vacua is diminishing and the depth of the former one is getting 
shallower. At h > 2 the electroweak minimum disappears at all. This pic- 
ture is typical and realized at any high temperature. Therefore, it is possible 
to obtain an upper limit on the magnetic field strength requiring that the 
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electroweak vacuum must be a long living state. From the above analysis 
one could conclude that the fields H > 2H = M 2 /e ~ 2 ■ 10 24 G had not 
been generated in the early universe. In the opposite case our world would 
never been realized and the system from the very beginning suppose to be in 
the global unbounded minimum. Similar symmetry behaviour (with slightly 
different values of h, B) has also been determined for K = 2. 

Fore completeness, let us describe the symmetry behaviour for the EP 
without the fermion contributions. In this case the symmetry restoration ( 
for K — 1,2) is also realized by the second type phase transition but typical 
temperatures to be of order B ~ 0.5 — 0.6. Naturally, no global minimum 
exists, so no limits on the magnetic field strength can be derived. 

To summarize the above results we stress that fermions affect in a very 
essential way the symmetry behaviour in the field H. We also recall that our 
consideration was based on the one-loop EP only. 

Now, let us include in our consideration the contribution of V r i ng . In Fig. 
6 we show the plot of V r i ng . It represents a complicated dependence on h 
and this contribution for strong fields acts to remove the separation barrier 
and stimulate the transition to the global vacuum state. 

In Figs. 7 ,8 and 9, 10 the influence of the ring diagrams is represented for 
small B , weak fields h ~ 0.01 — 0.1 and K = 0.5, 0.75 Respectively To better 
clarify their role we show the plots in parallel for chosen K. As is seen, with 
V T i n g included symmetry behaviour is considerably changed as compare with 
results presented in Fig. 5. Most important fact is that For K — 1 the local 
minimum is not realized at all even for weak field strengths. Actually, the 
value of K = 0.75 corresponding itlh ~ 70Gev is a bound value for the mass 
m H and the field strengths H ~ 0.01 - 0.1M 2 /e = 0.01 - 0.1 • 10 24 G give 
upper limits on the magnetic fields in the phase transition epoch. Stronger 
fields stimulates straight transition to the global minimum. Thus we see that 
due to V ring term the upper limit on the magnetic field is decreased from 
2 • 10 24 G derived with the EP V w to 0.1 • 10 24 G. Moreover, they restrict the 
Higgs boson mass: m# < 70Gev , otherwise the local electroweak minimum 
is not produced. We also observe that the phase transition is to be of the 
first type one. 

It is very important that the described minima of the EP are stable 
at high temperatures even for strong fields h. Really, typical value of B 
when the symmetry restoration happens for K ~ 1 is B ~ 0.1 . For the 
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calculated lower temperature mass we have j^g 2 / B 2 ~ 3. Hence we find that 
the effective mass of the unstable mode 2 — h + j^g 2 / B 2 is positive for all 
values of h considered. Thus, one has to conclude that classical constant 
magnetic field must be stable in the electroweak phase transition epoch. No 
W- and Z-boson condensates can be produced at high temperatures. These 
condensates would be realized at lower temperatures as the intermediate 
states of the vacuum. 



6. Discussion 

Let us discuss the vacuum stability in a magnetic field. This problem is 
of interest because of the presence in the jy-boson spectrum the mode 
pi = M 2 — eH which becomes unstable for H > M 2 /e. The evolution 
of this state and its consequences have been investigated in various aspects 



by many authors |25],[23|,[13]. At high temperatures, it was studied for the 
case when only the bosonic part of the electroweak theory has been included 
Efl,HI9|,|I5||. In Refs. [0 the classical equations of W, Z and fields were 



solved and the results that at high temperature the symmetry is restored and 
the magnetic field is screened by the W- boson and Z-boson condensates have 
been elaborated. The key point in this analysis is the assumption that sym- 
metry restoration is the second type phase transition which can be taken into 
account by adding the term ~ 2 T 2 in the field equations. For these results 
it was also important that the Higgs boson mass equals to rriz- From the 
point of view of the present investigation the described approach is not trusty 
because it does not take into account the vacuum polarization which is very 
important at high temperature and strong fields. In fact, this is the only one 
term produced by the vacuum polarization and other relevant contributions 
must be included. Exactly the polarization effects produce the temperature 
masses which stabilize the vacuum. So, no conditions for generation of the 
W- and Z-boson condensates exist. These condensates could be realized in 
strong magnetic fields at intermediate temperatures when imaginary part of 
the EP is nonzero. But in any case, to have a consistent picture the vacuum 
polarization should be taken into account because of a complicated behaviour 
of the EP in an intermediate range of T, H [|T2]],[iT?|. 



The present investigation carried out on the base of the complete EP 
with fermions included shown that at weak magnetic fields the symmetry 
restoration is to be the first type phase transition ( for the values of K 
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considered ,mn < M w ) . For strong fields our metastable vacuum could never 
be realized. But at any fixed temperature T there exists a corresponding 
magnetic field strength at which a metastable vacuum with positive energy 
and non-zero scalar field is produced . This picture may be of interest for 
cosmology. 

As follows from the results of Sect. 5 , the role of fermions is very essen- 
tial in the symmetry dynamics. Actually, their contribution determines the 
properties of the phase transition due to magnetic field. We have seen that 
at low temperature the field acts to prevent the phase transition from the 
metastable to stable minimum of the EP. Since the charged fermions and 
gauge bosons oppositely influence the symmetry behaviour in a magnetic 
field, for the actual values of particle masses these two contributions com- 
pensate each other and the metastable minimum position remains near the 
initial point <p 2 = 1 for any values of h. The influence of the field is expressed 
in the change of the potential barrier separating two minima. As is also oc- 
curred, at high temperatures the role of the ring diagrams is important (as 
also takes place at h — ||IS|,|I7[|). Thus, we conclude that the EW phase 



transition in a magnetic field acquires substantial changes. For its detailed 
investigation it is necessary to calculate the bubble nucleation parameters, 
the metastability bound on ms, etc. 

We would like to complete our discussion with a few remarks concerning 
the magnetic field in the restored phase. As was mentioned in the Intro- 
duction, in current literature a hypercharge magnetic field is discussed as a 
relevant one there [H] . There is an important difference between this field and 
ordinary magnetic field investigated in the present paper. The former field 
is an abelian one and requires an external source to maintain it in space. On 
the contrary, ordinary magnetic field may be considered as the projection of a 
nonabelian field produced spontaneously at high temperatures via Sawidy's 



mechanism. Investigation of the state was given in Refs.|7j|, p4[ ,||24j|. A final 
conclusion about this phenomenon can be done when the complete temper- 
ature mass of the unstable mode II (H, T) will be calculated. This work is in 
progress now. 
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Figure 3: Symmetry restoration at high temperatures and small magnetic 
fields determined by the one-loop effective potential V(0 2 , h, B, K). 
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Figure 4: Symmetry behaviour at zero temperature and strong magnetic 
fields determined by the one-loop effective potential V(0 2 , h, K). 
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Figure 5: Symmetry behaviour at fixed high temperature and a number of 
values h described by the one-loop efective potential V(</> 2 , h, B, K). 
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Figure 6: The V ring curves as the functions of (ft 2 for fixed temperature and 
a^nymber of h. 
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Figure 7: Symmetry behaviour at high temperature and 'weak' magnetic 
fields determined within the one-loop effective potential V = + V*- 1 - 1 for 
K = 0.5. 
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B=0.05 




Figure 8: Symmetry behaviour at high temperature and 'weak' magnetic 
fields determined within the total effective potential V = + + V r i ng 
%KT = 0.5. 
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Figure 9: Symmetry behaviour at high temperature and 'weak' magnetic 
fields determined within the total effective potential V = V^ -* + V^ 1 -* + V r i ng 
for K = 0.75. 



26 




27 



